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Abstract. I calculate the emission expected from a Poynting-flux-dominated gamma-ray burst (GRB) flow in 
which energy is dissipated gradually by magnetic reconnection. In this picture, the energy of the radiating 
particles is determined by heating and cooling balance (slow heating model). Detailed radiative transfer 
calculations show that, at Thomson optical depths of order of unity, the dominant radiative process is inverse 
Compton scattering. Synchrotron-self-absorbed emission and inverse Compton dominate in the Thomson thin 
parts of the flow. The electrons stay thermal throughout the dissipation region because of Coulomb collisions 
(Thomson thick part of the flow) and exchange of synchrotron photons (Thomson thin part). The resulting 
spectrum naturally explains the observed sub-MeV break of the GRB emission and the spectral slopes above 
and below the break. The model predicts that the 7-ray power-law tail has a high-energy cutoff typically in the 
~ 0.1 — 1 GeV energy range that should be observable with GLAST. The model also predicts a prompt emission 
component in the optical and UV associated with the GeV emission. Observations of the prompt emission of 
GRB 061121 that cover the energy range from the optical to ~ 1 MeV are explained by the model. 

Key words: Gamma rays: bursts - radiation mechanisms: general - methods: statistical 
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1. Introduction 

The GRB emission is the likely result of internal energy 
release in an ultrarelativistic flow. The dissipative and ra- 
diative mechanisms for the GRB largely remain uncertain. 
A popular model for the energy dissipation invokes inter- 
nal shocks in an unsteady flow (Paczynski & Xu 1994; 
Rees & Meszaros 1994). An alternative proposal is mag- 
netic dissipation in a strongly magnetized flow (Thompson 
1994; Spruit et al. 2001) 

Internal shocks are efficient in dissipating a large frac- 
tion of the kinetic energy of the flow provided that it is 
highly variable, i.e., composed of distinct ejection events 
with strong variation in their bulk Lorentz factor 7 (e.g. 
Kobayashi et al. 1997). Energy is dissipated by the shocks 
at the location of the collision of the shells. Particles are 
assumed to be accelerated on a very short timescale at the 
shock front to ultrarelativistic speeds and non-thermal dis- 
tributions. Subsequently, they radiate a fraction of the dis- 
sipated energy via synchrotron and inverse Compton pro- 
cesses. The relevant radiative mechanisms and the emitted 
spectra depend, to a large extent, on the shock micro- 
physics and the corresponding Thomson optical depth of 
the flow at the radius of the collision (see, e.g., Pe'er et al. 
2006) 

On the other hand, the energy dissipation that pow- 
ers the GRB emission may be gradual and distributed 
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over a large part of the volume of the flow. The en- 
ergy of the radiating particles is determined by the heat- 
ing/cooling equilibrium (slow heating model; Ghisellini & 
Celotti 1999; Stern & Poutanen 2004). Such an energy bal- 
ance is expected to lead to sub-relativistic or mildly rela- 
tivistic temperatures in the flow. Magnetic dissipation in 
a strongly magnetized flow can provide a possible physical 
setup where gradual dissipation is realized. As shown in 
Drenkhahn (2002) and Drenkhahn & Spruit (2002; here- 
after DS02), dissipation though reconnection takes place 
over several decades in radius; typically in both Thomson 
thick and thin conditions. Dissipation in the reconnection 
model is responsible for both the acceleration of the flow 
and the prompt emission. 



1.1. Emission from magnetic dissipation 

In previous works (Giannios 2006; Giannios & Spruit 2007; 
hereafter G06 and GS07 respectively), we studied the ra- 
diative transfer close to the Thomson photosphere of the 
flow in the context of the magnetic reconnection model. 
The detailed Monte Carlo calculations have shown that, 
due to energy release, the flow develops a hot photosphere 
with comoving electron temperatures of tens of keV. In the 
photospheric region. Coulomb collisions are fast enough to 
thermalize the electron distribution. Upscattering of pho- 
tons that are produced deeper in the flow by those hot 
electrons leads to a powerful photospheric emission; it ac- 



2 



Dimitrios Giannios: Prompt GRB emission from gradual energy dissipation 



counts for ~ 3 — 30% of the luminosity of the flow. The 
resulting E ■ f{E) spectrum has a characteristic ~ 1 MeV 
peak followed by a flat high-energy power-law emission. 

Here, I extend the radiative transfer calculation to the 
Thomson thin region of the flow where, in the reconnec- 
tion model, there can be substantial energy release and 
associated emission. Because of the strong magnetic fields 
synchrotron self absorption results in eSicient energy ex- 
change of the electrons, keeping their distribution thermal 
(Ghisellini et al. 1998). I show that heating and cooling 
balance in that region leads to electron temperatures of 
the order of the electron rest mass. Under these conditions, 
synchrotron-self-absorbed (SSA) emission is an important 
radiative mechanism in addition to inverse Compton. SSA 
dominates the observed emission in the soft X-rays and 
softer bands. 

The efficient thermalization of the emitting particles 

throughout the flow reduces the dependence of the model 
on the, poorly understood, mechanisms of particle accel- 
eration that operate in magnetic reconnection. The model 
is defined by just 3 main parameters (luminosity, baryon 
loading and a reconnection-rate parameter; see Sect. 3). In 
contrast to the internal shock model, no quantities have 
to added to parameterize the particle distributions and 
the amplification of magnetic fields; the field strength is 
an integral part of the reconnection model. Because the 
model is so well defined, it makes direct and stable pre- 
dictions for the emitted spectrum. In the following, I show 
how this spectrum can be computed. 

The structiirc of the paper is the following. In the next 
section I siimmarizc and contrast the main aspects of in- 
ternal shock and magnetic reconnection models for the 
prompt GRB emission. The dynamics of the fiow in the 
reconnection model and the radiative transfer in the flow 
are the topic of Sects. 3 and 4 respectively. The resulting 
spectra and their direct comparison with multi-frequency 
observations (from optical to 7-rays; see Page et al. 2007) 
of GRB 061121 are presented in Sect. 5. I discuss the re- 
sults and conclude in Sect. 6. 

2. Magnetic reconnection versus internal shocks 

The internal shock model for the GRB emission invokes 
high variability in the bulk Lorentz factor of the flow that 
leads to internal collisions. The location where the col- 
lision of two shells takes place depends on their initial 
separation and bulk Lorentz factors. These collisions can 
dissipate a substantial fraction of the kinetic energy of 
the flow. Internal shocks are assumed to lead to parti- 
cle acceleration and magnetic fleld amplification at the 
shock front. If electrons receive a large fraction of the 
dissipated energy then they are accelerated to ultrarela- 
tivistic speeds. They cool down radiatively by synchrotron 
and inverse Compton mechanisms. The resulting spec- 
trum and the relevant radiative mechanisms depend on 
details of the distribution of the accelerated particles, the 
magnetic field strength and the Thomson optical depth 
of the flow at the location of the collision. If the colli- 



sion of two shells takes place at the Thomson thin re- 
gion, as needed to explain the typical variability proper- 
ties of the GRB lightcurves (Daigne & Mochkovitch 1998; 
Nakar & Piran 2002; Mimica et al. 2005), synchrotron self 
Compton is likely the most promising radiative mecha- 
nism (Rces & Mcszaros 1994; Katz 1994; Tavani 1996). 
In this picture optically thin synchrotron emission domi- 
nates the observed hard X-ray, -^MeV spectrum. Despite 
its attractive features, the synchrotron model has theoret- 
ical and observational difficulties (as discussed, for exam- 
ple, in Ghirlanda et al. 2003). If, on the other hand, the 
collisions take place close to the Thomson photosphere, 
Compton scattering is the dominant radiative mechanism 
that shapes the spectrum and results in very different 
emission (Pe'er et al. 2006). 

As an alternative to internal shocks, magnetic dissi- 
pation can power the prompt emission provided that the 
flow is launched Poynting-flux dominated (or with a sub- 
stantial fraction of its energy in the form of Poynting 
flux). Magnetic dissipation through, for example, recon- 
nection can release energy smoothly in a large fraction 
of the volume of the flow. This energy release can take 
place while the flow expands over several decades in ra- 
dius (DS02). The energy of the radiating particles is deter- 
mined by balancing heating and radiative (or adiabatic) 
cooling at each radius. In this case the slow heating pic- 
ture described by Ghisellini & Celotti (1999) takes place. 
The electrons are subrelativistic or mildly relativistic and 
their synchrotron emission is self absorbed. In a strongly 
magnetized flow, such SSA emission guarantees efficient 
energy exchange and thermalization of the electrons on a 
very short timescale (Ghisellini et al. 1998). The resulting 
emission does not depend on details of particle acceler- 
ation and magnetic field amplification that one faces in 
internal shock models. The model is defined by just the 
luminosity, the baryon loading and the reconnection-rate 
parameter of the fiow. This contributes significantly to the 
predictive power of the magnetic reconnection model. The 
total observed fiux is the integrated emission from the dif- 
ferent parts of the flow in which dissipation of energy takes 
place. It contains both photospheric (Thompson 1994; 
Stern 1999; G06; GS07) and Thomson thin components 
(this work). 

2.1. Implications from the observed GRB variablity 

One additional difference of the internal shock and the 
magnetic reconnection model is connected to implications 
from the observed variability of the lightcurves. Internal 
shocks are efficient only in variable flows. Variability and 
dissipation are, a priori, unrelated in the magnetic re- 
connection model in which dissipation takes place, even 
in a steady outflow. On the other hand, the observed 
lightcurves are often highly variable showing that the flow 
does evidently evolve during a GRB. In the context of 
the reconnection model, the observed variability reflects 
changes in the luminosity and baryon loading of the flow 
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during the burst. As shown in GS07, the flow can be 
treated as quasi-stationary for all but the shortest time 
scales observed in a burst, the variation of spectral proper- 
ties during a burst directly reflects variations in the central 
engine. This is in contrast to models in which the prompt 
radiation is produced at much larger distances from the 
source, such as external shock models. It is also in contrast 
with the internal shock model, since the internal evolution 
of the flow between the source and the level where radia- 
tion is emitted is a key ingredient in this model. Deducing 
properties of the central engine from observed burst prop- 
erties is thus a much more direct prospect in the magnetic 
dissipation model. 

3. Gradual energy release because of magnetic 
reconnection 

Magnetic dissipation can take place gradually in the GRB 
flow. The rate at which energy is dissipated as a function of 
radius depends on the magnetic field geometry and the ex- 
act mechanism through which magnetic energy dissipates. 
If the flow is launched with field of large scale, energy dis- 
sipation can be a result of global MHD instabilities (such 
as current-driven instabilities; e.g., Lyutikov & Blandford 
2003; Giannios & Spruit 2006). On the other hand, if the 
flow contains reversing magnetic fields of sufficiently small 
scale, dissipation can take place directly through recon- 
nection (Drenkhahn 2002; DS02; Thompson 2006). Here, 
I focus on the reconnection model, which makes clear pre- 
diction for the energy dissipation as function of radius; 
essential for the radiative transfer calculations presented 
here. Though the results presented here are directly appli- 
cable to the DS02 model, qualitatively similar results are 
expected from other gradual, magnetic dissipation models. 

3.1. The reconnection model 

An important physical quantity in the reconnection model 
is the ratio ctq of the Poynting flux to kinetic energy flux 
at the Alfven radius tq. This quantity parameterizes the 
baryon loading of the flow j] and determines the terminal 

3/2 

bulk Lorentz factor of the flow ~ r/ ~ Cg . The flow 
must start Poynting-flux dominated with cto ^ 30 for it to 
be accelerated to ultrarelativistic speeds with 700 100 
that are relevant for GRB flows. 

In the reconnection model, the magnetic field in the 
flow changes polarity on a scale A. If the magnetic field 
anchored in the rotating central engine is nonaxisymmet- 
ric, this scale is (in the central engine frame) of the order 
of the light cylinder n: A ~ 27rc/il, where is the angular 
frequency of the rotator. This is as in the oblique rotator 
model for pulsar winds (Coroniti 1990). This model has 
been further developed to include special relativistic ef- 
fects and different reconnection prescriptions (Lyubarsky 
& Kirk 2001; DS02; Kirk & Skjaeraasen 2003). The rate 
of magnetic reconnection DS02 model is parameterized 
through the velocity with which magnetic fields of op- 
posite direction merge. The Vr is assumed to scale with the 



Alfven speed, wa, i.e. Vr = evA- A nominal value used for e 
is 0.1 (see Lyubarsky 2005). For the flows with (To ^ 1 that 
are of interest here, the energy density of the magnetic 
field is larger than the rest mass energy density, hence 
VA ~ c, and the reconnection takes place with subrela- 
tivistic speeds. 

3.1.1. Properties of the flow 

In the reconnection model, magnetic dissipation takes 
place all the way from the initial radius tq till the sat- 
uration radius r^. Part of the dissipated energy (approx- 
imately half) is directly used to accelerate the flow. The 
acceleration of the flow is gradual following the 7 ~ r^/"^ 
scaling as function of radius in the regime tq <S r <^ Ts. 
To first order approximation, no further acceleration takes 
place beyond the saturation radius. Summarizing the re- 
sults derived in Drenkhahn (2002), the bulk Lorentz factor 
of the flow is approximately given by 
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7 = 7oo 1 - ) =148 rli^ien)l^^al/l for r < 



(1) 



7 = 7oo = o■o^^ for r > r^. 
while the saturation radius is 



7rc7go 
Sen 



or rs.ii = 310 



"0,2 



(2) 



The notation A — 10^ is used; the 'reference values' 
of the model parameters are ctq = 100, e = 0.1, 17 = 10^ 
rad-s~^. The product of e and fl parameterizes the recon- 
nection rate. The physical quantities of the flow depend 
on this product. 

In the steady, spherical flow under consideration the 
comoving number density can be written as 



L 
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(3) 
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where L is the luminosity per steradian of the GRB flow. 
The reference value used is L = 10^^ erg-s~"'^-sterad~^. (In 
this form the expression can be compared with the fireball 
model, where the baryon loading parameter 77 replaces the 
factor ctq^^). 

The expression ([T]) is deviating from the exact numeri- 
cal solution presented in Drenkhahn (2002) at r > r^. The 
reason is that the dissipation does not stop abruptly at 
Tg but there is modest energy release at a slower rate at 
larger radii. This leads to modest acceleration of the flow 
at r > Tg. In the following, we ignore these deviations and 
use the expressions (1) and (3) for the bulk Lorentz factor 
and density of the flow respectively. This simplification 
facilitates the radiative transfer study in the flow. 

On the other hand, quantities of the flow such as mag- 
netic field strength and the energy dissipation rate as func- 
tions of radius need to be followed in more detail around 
Tg. Though, not important for the global energetics, the 
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remaining dissipation at the radii r > Tg results in syn- 
chrotron emission that dominates the observed radiation 
in soft bands (such as optical and near ultra violet). This 
emission is mainly a result of the large emitting surface 
at these outer parts of the flow. I take into account the 
residual dissipation, so that to correctly describe the soft 
emission. In the calculations that follow I calculate the 
magnetic field strength and dissipation rate by solving 
the relevant differential equations that describe the full 
ID relativistic MHD problem as presented in Drenkhahn 
(2002; summarized by Eq. (38) in that paper). Still, for the 
purpose of estimates, I give analytic expressions for these 
quantities that are accurate below the saturation radius. 

The comoving magnetic field strength B' below the 
saturation radius is given by setting L ~ Lp = 
r^B'^j'^c/i-rr and solving for B': 



B' = 



AttL \ 1/2 



9 9 
CT I 



for r < i\ 



(4) 



The rate of energy density release in a comoving frame can 
be found by the following considerations. The time scale 
over which the magnetic field decays is that of advection 
of magnetic field of opposite polarity to the reconnection 
area. The reconnection speed is Vr = sva — £c, while the 
magnetic field changes polarity over a comoving length 
scale A' = 2Trjc/n. The decay timescale for the magnetic 
field, therefore, is 



_ A' _ 27r7 

^dcc 7Z~ ■ 

Vr Sll 



(5) 



Using the last expression and Eqs. ([T]) and the rate 
of dissipation of magnetic energy density in the comoving 
frame is 



-^diss 



(Br 



em 



idcc/2 27rcr273 



for r < r. 



(6) 



The bulk Lorentz factor 7, the density n', the magnetic 
field strength B' and rate of dissipation of magnetic energy 
density Pdiss of the flow as functions of radius are the 
quantities needed for the study of the resulting emission. 

4. Photospheric and Thomson thin emission 

If all the energy is dissipated deeply into the GRB flow 
(i.e. at large optical depths), adiabatic expansion converts 
most of this energy into kinetic at the expense of radia- 
tion. Gradual dissipation heats the flow continuously and 
maintains a substantial fraction of the energy in the form 
of radiation. This radiation is released at the photosphere 
of the flow. If dissipation takes place further out in the 
flow it can result in additional emission coming from the 
Thomson thin region. The total flux received by the ob- 
server is the integrated emission from the different parts 
of the flow where dissipation takes place. 

In the case of magnetic dissipation (as well as for other 
dissipative mechanisms), the fate of the released energy is 
rather uncertain. An interesting possibility is that dissipa- 
tion leads to MHD turbulence where particle acceleration 



can take place by scattering of photons by Alfven waves 
(Thompson 1994). On the other hand, the magnetic en- 
ergy can directly be dissipated to the particles in the flow, 
most likely to the electrons due to their higher mobility. 
Following G06, GS07 we assume that a fraction /o of order 
of unity of the dissipated energy heats up the electrons. 
For the results presented here we set /e = 0.5. 

The resulting emission does not depend only on the 
amounts of energy released but also on the distribution of 
the emitting particles. I assume that the electron distri- 
bution is thermal throughout the region where dissipation 
takes place. As I discuss in more detail in Sect. 4.2.1., 
the thermalization of the electrons is result of Coulomb 
collisions in the inner parts of the flow and of exchange 
of synchrotron photons at the outer parts. I first summa- 
rize the results of G06, GS07 on the photospheric emission 
from the reconnection model and then turn to the study 
of the Thomson thin emission. 



4.1. Photospheric emission 

In addition to the saturation radius rg, another charac- 
teristic radius of the flow is the Thomson photosphere. 
The Thomson optical depth as a function of radius is 
T ~ n'aTr/j. It can be expressed in terms of the pa- 
rameters of the fiow (e.g. G06): 
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L 



52 
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(7) 
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As expected, at small radii to optical depth is large and 
vice-versa. The radius of the Thomson thick-thin transi- 
tion is found by setting t = 1 in Eq. (7) and solving for 
rph- 



rph.ii 



= 6 



L 



3/5 

52 



/ o\2/5 3/2" 

(en)^' <j' 



(8) 



In deriving these expressions, we have assumed that rph < 
Tg. A similar calculation gives the radius of the photo- 
sphere in the rph > case. 

One can check that for a large parameter space relevant 
for GRB flows, rph < which means that dissipation 
proceeds throughout the photospheric region. In terms of 
the physical properties of the flow, there is a critical value 
of the magnetization cro,cr, f^or which rph — rg. For ctq > 
<^o,cT, Tph < rs- Using Eqs. ^ and ^ one finds 



cro,cr = 42(i52(el7); 



V 2/15 



(9) 



The critical baryon loading depends weakly on the param- 
eters of the flow: rjcr — ctq „ = 270(L52(eri)3)^^^. 

For (To ^ fccr dissipation ceases deep in the flow (at 
high optical depths). For (Tq S> fo^cr, dissipation takes 
place in both Thomson thick and thin conditions with 
most of the energy released in the outer parts of the flow. 
The implications of such energy release to the properties 
of the flow and the resulting radiation have been studied 
in G06 and GS07. Those studies focused on the Thomson 
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thick part of the flow and the photospheric region out to 
T ~ 0.1. Observational effects from Thomson thin dissi- 
pation were not considered in detail and are the topic of 
this paper. 

The main results of the GOG, GS07 papers are the fol- 
lowing. Particles and radiation are found to be in ther- 
mal equilibrium deep in the flow. There, the comoving 
temperature Tth of the flow is calculated, under the as- 
sumption of complete thermalization, by integrating the 
energy released at different radii in the flow and taking 
into account adiabatic cooling. Due to the dominance of 
scattering, the details of radiative transfer become im- 
portant already at fairly large optical depth in the flow. 
Equilibrium between radiation and matter holds only at 
Thomson depths greater than about 50. At smaller op- 
tical depths the electron distribution stays thermalized, 
but is out of equilibrium with the photon field. More dis- 
cussion on the processes that lead to thermalization of 
the electron distribution is presented in the next section. 
Compton scattering of the photons is treated in detail in 
this region with Monte Carlo Comptonization simulations 
(G06; GS07). Energy dissipation at moderate and low op- 
tical depths is shown to lead to emission that has a highly 
non-thermal appearance. 

For (To ~cro,cr, the flow develops a hot photospheric re- 
gion. The electron temperature at moderate optical depths 
can be estimated analytically by balancing the heating 
[Eq. ^] with the inverse Compton coohng rate Pcomp = 
AkBTcu'caTUi/nicC^, where u^- is the energy density of ra- 
diation. The electron temperature can be expressed as a 
function of optical depth in the flow (the detailed numer- 
ical results verify these estimates; see G06): 

To~— keV, for 0.1<r<50. (10) 

T 

For T > 50, the electron temperature equals that of the 
photon field. For r < 0.1, the electrons become mildly rel- 
ativistic and one has to consider relativistic corrections 
to the Compton cooling to derive the appropriate expres- 
sion for the temperature. At these larger radii (or lower r) 
the electron temperature becomes high enough that syn- 
chrotron emission (and associated cooling of the particles) 
has to be taken into account. 

4.2. Emission from the Thomson thin region 

Magnetic dissipation around the Thomson photosphere 
leads to subrelativistic electrons that are more energetic 
with respect to the average photon in the flow. Under these 
conditions the dominant mechanism that shapes the spec- 
trum is inverse Compton scattering. Synchrotron emission 
is negligible since it is strongly self absorbed. When dis- 
sipation continues at the optically thin parts of the flow 
the electron temperature becomes of the order of the elec- 
tron rest mass. At those temperatures, synchrotron self 
Compton emission has important effect on the emitted 
spectrum (see also Stern & Poutanen 2004). Furthermore, 
SSA affects the electron distribution in the ffow. 



4.2.1. Thermalization of the electrons 

At high optical depths, the density of the ffow is high 
and the electron temperature rather low [see Eq. (fTO|) ]. 
Under these conditions. Coulomb collisions are efficient 
in thermalizing the electron distribution. This is not the 
case in the Thomson thin parts of the flow (G06). On the 
other hand, there is a more efficient channel for energy 
exchange among electrons. It is shown in Ghisellini et al. 
(1998) that electrons with energy of the order of their rest 
mass can thermalize on a few synchrotron cooling times 
(as defined for thin synchrotron emission) by emitting and 
absorbing synchrotron photons. In the magnetically dom- 
inated flow under consideration, the thermalization of the 
electrons takes place on a timescale shorter than the heat- 
ing/cooling one. One can, therefore, assume that the elec- 
trons are approximately thermal when calculating their 
emission. 

A caveat in the previous argument is connected to the 
underlying assumption that the dissipated energy is dis- 
tributed to a large fraction of the particles in the flow. If 
the magnetic reconnection leads, for example, to deposi- 
tion of most of the energy to a small fraction of the elec- 
trons, they can be accelerated to relativistic speeds and 
cool efficiently through thin synchrotron emission. In this 
case, synchrotron self absorption is not an efficient ther- 
malization mechanism. This case has been investigated 
in Giannios & Spruit (2005). From this point on we as- 
sume that the dissipated energy is distributed among a 
large fraction of the particles and hence thermalization is 
achieved in the electron distribution. 

4.2.2. Modeling of the synchrotron emission 

The electron temperature becomes mildly relativistic at 
small optical depths r ~ 0.1 (see Eq. ([TU])'). It increases 
further at larger radii resulting in substantial synchrotron 
emission. The synchrotron-self-absorbed emission from 
mildly relativistic plasma has a characteristic spectrum 
that consists by a Rayleigh- Jeans part up to the so-called 
turnover frequency vt where the optical depth due to syn- 
chrotron absorption of the flow becomes unity. Most of 
the energy is emitted at the turnover frequency. At higher 
frequencies, the spectrum is very steep following the ex- 
ponentially decaying tail of the synchrotron thin emis- 
sion. The energy density per unit frequency of the syn- 
chrotron photons in a frame comoving with the flow is 
— Snv^kBT/c^ for u < vt, while it drops very fast for 
v > Vf The energy density of the synchrotron photons is 
given by integrating the last expression: Us = J Uj^dv ~ 
STTv^ksT/'ic'. The synchrotron luminosity per steradian 
at radius r is: 

The synchrotron emission depends strongly on the 
turnover frequency. The turnover frequency can be related 
to the magnetic field strength B' , the electron temperature 
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Te and the Thomson optical depth r of the flow. For the 
mildly relativistic plasma under consideration the calcu- 
lation of the synchrotron emission and the corresponding 
absorption is rather evolved. On the other hand, there are 
studies focusing on developing approximate expressions 
for the synchrotron emission (and therefore absorption) at 
this temperature regime (e.g. Petrosian 1981, Mahadevan 
& Narayan 1996; Zdziarski et al. 1998; Wardzihski & 
Zdziarski 2000). Here I use the approximate analytic ex- 
pressions for i>t derived in Zdziarski et al. (1998): 

^^t, = ^0e'^cln^^4-, (12) 



36 



In 



(^c ~ k-QTc/iricC^. The quantity C is defined a^j 



where — eB' jlixmcC is the cyclotron frequency and 



(13) 



C 



3 

7^ 



2/7 



where h, af are Planck's and fine structure constants re- 
spectively. For typical model parameters and radii r ^ Vg, 
the turnover frequency is Vt ^ lO'^J^c which results in 
z/°bs =7^i/j/(l + z)<l keV {z is the redshift of the burst). 
The synchrotron-self-absorbed emission appears mainly in 
the soft X-rays and softer bands. 

The expression ([T2|) has been compared with the ex- 
act numerical results and shown to be overestimating 
the turnover frequency (and consequently the synchrotron 
emission) in plasma with 0c ^0.1 (see fig. 2b in Wardzihski 
& Zdziarski 2000). On the other hand, it is quite accu- 
rate for higher electron temperatures. Since most of the 
synchrotron emission comes from regions of the flow with 

~ 1, Eq. (|12p is accurate enough for the calculations 
presented here. 

One can now extend the Comptonization calculation 
developed in G06 to include the radiative transfer in the 
Thomson thin region of the flow. The procedure is the 
following. I make a choice for the temperature profllc as 
function of radius of the electron temperature in the flow. 
Analytical expressions such as that of Eq. (10) provide 
a good initial guess. The radiative transfer in the flow 
is studied using the Monte Carlo Comptonization code 
described in G06. The calculation includes the thermal 
radiation field carried with the fiow which is injected in 
the inner numerical boundary at the "equilibrium" radius 
where radiation and particles drop out of equilibrium. At 
larger radii the synchrotron emitted flux is also included 
(using expressions (fTTj) . (fT2|) ). Both sources of photons are 
propagated through the medium and their scattering by 
electrons is followed. The code calculates the spectrum 
and the radius-dependent cooling rate of the electrons. 
Cooling because of inverse Compton, synchrotron emission 
and adiabatic expansion is taken into account. The adi- 
abatic cooling becomes important at high optical depths 
and at the very outer parts of the fiow where the expan- 
sion timescale r/jc is shorter than the radiative cooling 



^ Here I set the Zdziarski et al. (1998) correction factor A-\ 



one. The outer boundary of the calculations is set at large 
enough radius so that it does not have an effect on the 
computed spectra. I iterate the electron temperature un- 
til the cooling rate matches the heating rate predicted by 
the model reasonably well at all radii. In practice, I make 
sure that they match within ^ 30% or less everywhere. 
Considering the rather larger uncertainties in the model 
coming from, for example, the assumed dynamics of the 
fiow, this a fairly accurate calculation. 

5. Resulting spectra 

First, I focus on the new features that appear in the emit- 
ted spectrum w.r.t. the G06 results because of the inclu- 
sion of the Thomson thin emission. In the illustrative case 
of Fig. 1, I set (To = 70 (which corresponds to fiow with 
baryon loading rj ~ CTq^'^ 600) and the rest of the pa- 
rameters to their reference values. Spectra are plotted in 
the central engine frame. The dotted line shows the input 
thermal radiation at the "equilibrium radius" which is the 
inner boundary of the computed domain. The thermal fiux 
is advected with the flow and constitutes a large fraction 
of the seed photons to be Comptonized further out in the 
fiow. 

The appearance of the photon spectrum at radius 
which corresponds to optical depth r = 0.1 is shown with 
dashed line. This radius is the outer boundary used in 
most of the calculations of G06, GS07. One can clearly 
see the effect of inverse Compton scattering to the spec- 
trum. The peak of the E- f{E) spectrum increases slightly 
and the spectrum becomes broader. An important fea- 
ture is the high-energy tail that is result of unsaturated 
Comptonization at r < 1. Note also a second peak at 
'-^l kcV. This is result of synchrotron emission with the 
turnover frequency being '^l keV (in the central engine 
frame). This component is still weak at r = 0.1 and has 
not been included in the calculations of G06, GS07. 

The total spectrum is shown with solid line. This in- 
cludes the emission from the whole volume of the flow 
where there is dissipation taking place. The overall emis- 
sion spectrum is much broader. It is characterized by 
a break at ~ 1 MeV followed by a fiat spectrum with 
photon-number index F ~ —2 (where dN/dE ~ E^), 
close to the typically observed one. The hard 7-ray tail 
is extending up to ~1 GeV which corresponds to the 
Lorentz boosted temperature of the flow at its outer lay- 
ers (where it reaches its maximum values) . Comparing the 
spectrum at r = 0.1 and the total one, it is clear that 
Comptonization proceeds throughout the Thomson thin 
region strengthening the hard 7-ray component. 

An important new feature is the powerful emission 
that appears in the soft X-rays and softer bands. This 
comes from synchrotron-self-absorbed emission. SSA dom- 
inates by many orders of magnitude the ultra violet and 
optical emission. The softer emission originates from the 
Thomson thin part of the fiow that is characterized by 
the higher electron temperatures and larger emitting sur- 
face (see also Stern & Poutanen 2004). This emission is 
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Fig. 1. Photon spectrum at different radii (or correspond- 
ing Thomson optical depths) in the flow. The spectrum is 
shown in the central engine frame. The dotted line stands 
for the spectrum at the radius where radiation and elec- 
trons decouple. The photospheric emission is shown with 
the dashed line (see also GOG; GS07). The overall spec- 
trum (solid line) includes the emission from the Thomson 
thin region of the flow. Synchrotron-self-absorbed emis- 
sion from this region dominates the spectrum below ~ 10 
keV. Inverse Compton leads to flat 7-ray tail up to ~ 1 
GeV. 

very weak in models where dissipation takes place be- 
low or around the Thomson photospheric region (see, e.g., 
Meszaros & Rees 2000; Pe'er et al. 2006; loka et al. 2007). 

As a result of the synchrotron-self-Compton compo- 
nent, the spectrum below the MeV peak softens with re- 
spect to the G06 calculation where only the photospheric 
component was considered. The spectrum can be well fit- 
ted with a power-law in the 30 — 300 keV energy range 
with photon- number index of F ~ —1.2 which very close 
to the one typically observed (e.g. Preece et al. 1998). 

In this example, the SSA emission spectrum hard- 
ens considerably below ^ 30 eV. At lower energies the 
Rayleigh- Jeans limit is gradually approached. The loca- 
tion of the hardening is determined by the radius where 
adiabatic expansion starts to dominate the cooling of the 
electrons. In this example adiabatic expansion dominates 
at r > Tad ~ 5 • 10^^ cm. Most of the radiation observed 
at E < 30 eV comes from this radius. The optical and 
near UV emission is delayed w.r.t. the ~ MeV emission by 
St ~ TsAllto^ ~ 0-05 s. Radiation in these bands reaches 
the observer with small but maybe detectable lags. In the 
region 30 < i? < 1000 eV the spectral slope depends on the 
radial dependence of the turnover frequency v\_. 

The relative strength of the synchrotron-self-Compton 
(SSC) component depends on the fraction of energy dis- 
sipated in the Thomson thin region of the flow. For 
(Jo ^ cro,cr most of the energy is dissipated in the Thomson 
thin region. Correspondingly the SSC component is pro- 
nounced. This is evident in Fig. 2 where the spectrum is 



Fig. 2. Resulting spectrum (in the central engine frame) 
for different baryon loadings of the flow. From bottom to 
top the curves correspond to magnetization cto =40, 50, 
60, 70, 100 (or corresponding baryon loading t] ~ 250, 
350, 460, 590, 1000) respectively. The high cto flows are 
characterized broader spectra. The model predicts that 
bright prompt optical and UV emission is accompanied 
by powerful ~GeV emission. For bright optical emission, 
the optical spectrum is expected to be hard. 



shown for different values of cto (the rest of the parameters 
are kept in their reference values). For ctq = 40, dissipa- 
tion stops close to the photosphere of the flow. The SSA 
component is almost absent and the emission above the 
thermal peak at ~ 1 MeV relatively weak. In more baryon 
loaded models the emission is quasi-thermal since dissipa- 
tion stops at high Thomson depths. 

With increasing ctq both the SSA and inverse Compton 
components become relatively more powerful. The ther- 
mal peak is followed by a flat 7-ray emission that extends 
up to ~0.1-1 GeV. Most of the models show a high energy 
cutoff in this energy range. This cutoff corresponds to the 
highest energies to which photons are upscattered. It is 
determined by the Lorcntz boosted electron temperature 
at r ~ T-g. The spectral slope below the ~1 MeV break 
becomes softer with increasing ctq. The photon- number 
index in the 30 — 300 keV energy band varies in the range 
— 1.2 < F < —0.4 in agreement to that typically observed 
(e.g. Preece et al. 1998). The high ctq models have power- 
ful optical and near ultra violet emission. The flux /(i?) 
that is emitted in these bands is similar to the X-rays one. 
The optical spectrum is hard with photon number index 
<F< -f-1. 

Varying the baryon loading of the flow has moderate 
effect in the emission in the BA TSE energy range but pro- 
foimd implications in other bands. The model predicts 
that flows with low baryon loading (i.e. high ctq) have 
powerful optical, UV and GeV emission. More on the com- 
parison of the model with observations is presented in the 
next section. 
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5.1. Comparison with observations 

The prompt GRB emission has been typically observed in 
the hard X-rays up to ~ 1 MeV 7-rays. The spectrum in 
this energy range shows a characteristic sub-MeV break 
followed by a flat power-law 7-ray tail (e.g. Band et al. 
1993). Below the break the spectrum has typical photon 
number index of F ^ —1 although much harder spectra 
have also been observecO. The observed sub-MeV break 
and the spectral slopes above and below the break are 
naturally explained by the gradual dissipation model dis- 
cussed here. Furthermore, the model makes specific pre- 
dictions on the prompt emission from the optical to GeV; 
bands that are currently (or will soon be) accessible to 
observations. The model predicts that the flat 7-ray tail 
extends up to a cutoff that typically appears at ~ 0.1 — 1 
GeV. It also predicts the prompt optical and UV emis- 
sion. For low baryon loading, the emission in these bands 
is powerful with energy flux f{E) similar to that of X-rays. 
The optical emission is characterized by a hard spectrum. 
Optically bright bursts have powerful GeV emission and 
softer spectra below the ~1 MeV break. 

In recent years, several observations in softer bands 
have been made simultaneously with the prompt GRB 
emission. The Swift satellite has observed the prompt 
emission in the X-rays (Hill et al. 2006; Romano et al. 
2006; Page et al. 2007) and uhra violet (Page et al. 2007) 
and robotic telescopes in the optical and infra red (e.g. 
Akerlof 1999; Vestrad 2005; Blake 2005; Boer et al. 2006; 
Vestrand et al. 2006; Klotz et al. 2006). Furthermore, 
GLAST is expected to probe the emission from GRBs up 
to - 100 GeV. 

Here, I compare the model to the very well sampled 
prompt emission of GRB 061121. This burst has been ob- 
served from optical to ~ 1 MeV (Page et al. 2007). The 
prompt emission has been followed with XRT and UVOT 
on board to Swift and in the optical with iJOTS'i? simulta- 
neously to 7-ray observations with BAT and Konus- Wind. 
There are two time resolved spectra just before and dur- 
ing the main pulse of the prompt emission that appears 
~ 75 sec after the onset of the burst. The pulse is clear in 
the lightcurves in all observed energy bands. The correla- 
tion between the different bands indicates that the optical, 
UV and X-ray and 7-ray components have common ori- 
gin (i.e. they are connected to the prompt GRB). This is 
unlike cases where the optical lightcurves are not tracing 
the 7-rays (e.g. Akerlof 1999; Boer et al. 2006; Klotz et al. 
2006 ) suggestive of a different physical origin w.r.t. that 
of the prompt emission. 

In Fig. 3, the data shown with circles refer to the pre- 
spike emission (epoch I in the Page et al. 2007 terminol- 
ogy) and the stars to the peak observed luminosity of the 
burst (epoch II). The data span approximately 6 orders 
of magnitude in frequency from the optical to 1 MeV. 
Overplotted are the spectral predictions of the model for 
two different sets of parameters. The low luminosity model 

^ These hard spectra cannot be explained by the thin syn- 
chrotron model (e.g. Ghirlanda, Celotti & Ghisellini 2003). 
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Fig. 3. Applying the model to multi frequency observa- 
tions of the prompt emission of GRB 061121 (see Fig. 11 
in Page et al. 2007). The circles stand for the observations 
of epoch I (just before the main pulse of the burst) and the 
stars for those of epoch II (during the pulse) . Observations 
are blue-shifted by 1+z to the burst rest frame {z = 1.131 
for GRB061121). The solid and dashed curves show spec- 
tra for two different sets of the parameters of the flow that 
illustrate that the model can account for the broad-band 
prompt spectra. 

has L = 10^° erg-s~^-sterad~^ and do — 38 and the high 
luminosity one L = 10^^ erg-s~^-sterad~^ and uo = 60. 
The two models (not meant to be detailed fits) are repro- 
ducing the observations quite closely. 

Note that for a given observed luminosity, the baryon 
loading is essentially the only free parameter of the model. 
This can be constrained by the ratio of the ^ IMeV-to- 
optical flux. Additional constraints can come from obser- 
vations of the prompt emission in harder bands. In this 
respect GLAST observations in the ~GeV range are go- 
ing to be of particular importance. 

The high luminosity model (that describes the epoch II 
observations) is characterized by higher (Tq with respect to 
the lower luminosity one. This is in qualitative agreement 
with the baryon loading-luminosity correlation during the 
evolution of the burst needed to explain observed energy- 
dependent properties of the GRB pulses in the context of 
the reconnection model (for details see sect. 4 in GS07). 
However since GS07 do not consider the Thomson thin 
emission in the calculations, the quantitative results of 
section 4 in GS07 have to be revisited. 

6. Discussion/conclusions 

The GRB emission may be result of internal collisions in 
a variable flow (Rees & Meszaros 1994). In the internal 
shock model energy is dissipated at the location of the 
shell collision. Particles are accelerated at the shock front 
to ultrarelativistic speeds and non-thermal distributions. 
Magnetic fields are assumed to be amplified because of 
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plasma instabilities. The emitted spectrum depends on the 
distribution of the emitting particles and the strength of 
the shock amplified magnetic field, both of which are not 
understood from first principles. The emitted spectrum 
also depends on the properties of the flow (such as optical 
depth) at the radius of the collision (e.g. Pe'er et al. 2006). 

As an alternative to internal shocks, magnetic dis- 
sipation in a strongly magnetized flow can power the 
GRB (Thompson 1994). Magnetic dissipation may lead 
to gradual release of energy over a wide range of radii 
(e.g. Drenkhahn 2002; DS02). It typically proceeds in both 
Thomson thick and thin regions of the flow. The released 
energy can be distributed to a large fraction of the parti- 
cles of the flow leading to the slow heating scenario for the 
GRB emission (Ghisellini & Celotti 1999). The emitting 
particles (i.e. electrons) are heated up to mildly relativistic 
speeds. Because of the strong magnetic fields exchange of 
synchrotron photons provides an efficient mechanism for 
the thermalization of the electron distribution (Ghisellini 
et al. 1998). Since the emitting particles are thermal the 
resulting emission does not depend sensitively on poorly 
understood physics of particle acceleration in magnetic re- 
connection. The model is defined by just the luminosity 
of the flow, its baryon loading and a reconncction-rate pa- 
rameter and makes direct and stable predictions for the 
electromagnetic spectrum. 

In previous works (G06 and GS07), we calculated 
the photospheric emission expected from the reconnection 
model. The radiative transfer study was made with Monte 
Carlo simulations. Those calculations have shown that the 
flow is characterized by powerful photospheric emission 
with most of the energy appearing in the hard X-rays and 
~ 1 MeV 7— rays. This emission is the result of photons, 
produced deep into the flow, that are inverse Compton 
scattered by sub-relativistic electrons at Thomson optical 
depths of order of unity. 

In the reconnection model, for low enough baryon- 
loading of the flow, there are substantial amounts of en- 
ergy dissipated in its outer, Thomson thin, parts. Here 
I have extended the G06 calculations to include the 
Thomson thin emission. Energy release at large radii leads 
to mildly relativistic electrons that cool down though 
emitting synchrotron radiation and inverse Compton scat- 
tering soft photons. SSA emission dominates the observed 
radiation in the soft X-ray and softer bands. This soft (e.g. 
optical) emission is very weak in models where dissipation 
is limited below or around the Thomson photosphere (see, 
e.g., Meszaros & Rees 2000; Pe'er et al. 2006; loka et al. 
2007). Inverse Compton in the Thomson thin region leads 
to a flat high-energy spectrum that extends up to GeV 
energies. 

The resulting spectra from the radiative transfer cal- 
culations naturally explain the observed sub-MeV break 
of the GRB emission and the spectral slopes above and be- 
low the break (Band et al. 1993). Furthermore, the model 
makes rather robust predictions for the emission in other 
energy bands. The flat 7-ray spectrum is expected to show 
a cutoff in the ~ 0.1 — 1 GeV energy range that should 



be observable with GLAST. The optical and ultra violet 
emission can be powerful and the optical spectrum hard 
with photon number index < F < 1. Bright prompt opti- 
cal emission is predicted to be accompanied by powerful ~ 
GeV emission and rather soft spectra below the sub-MeV 
break. Comparison with multi-frequency observations of 
the the prompt emission from GRB 061121 that span from 
the optical to the ~ MeV range (Page et al. 2007) supports 
the model. 

Acknowledgements. I thank Henk Spruit for many valuable 
suggestions and discussions. 
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